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Atmospheric neutrino-fluxes are calculated over the wide energy range from 30 MeV to 3,000 GeV 
for the study of neutrino-physics using the data from underground neutrino-detectors. In this 
calculation, a full Monte Carlo method is employed for low energy neutrinos (30 MeV - 3 GeV), 
while a hybrid method is used for high energy neutrinos (1 - 3,000 GeV). At low energies, the ratio 
(y e + v e )/(v fi +Vfi) agrees well with other calculations and is significantly different from observations. 
For the high energy neutrino-fluxes, the zenith angle dependence of atmospheric neutrino-flux is 
studied in detail, so that neutrino-oscillation parameters can be calculated for comparison with 
experimental results. The atmospheric muon-flux at high altitude and at sea level is studied to 
calibrate the neutrino-fluxes at low energies and high energies respectively. The agreement of our 
calculation with observations is satisfactory. The uncertainty of atmospheric neutrino-fluxes is also 
studied. 



I. INTRODUCTION 



> 

On 

m 
m . 

In this paper, we report the calculation of the atmospheric i/-flux in the energy range from 30 MeV up to 3,000 GeV, 
corresponding to the observation range of underground neutrino detectors. Detailed calculations of atmospheric in- 
fluxes are important, since the flux ratio, u^/y e , observed by many experiments shows a significant deviation from 
' the expected value |S| at low energies (< E v >^ 1 GeV) and multi-GeV energies (< E v >~ 5 — 7 GeV) [|. 
Many authors have considered the possibility that this deviation is evidence for ^-oscillations, with a large mixing 
O ' angle and Am 2 ~ 10~ 2 eV 2 [jjj |] The zenith angle variation of the (/Ve)data/(/V e )MC ratio at multi-GeV 
energies is especially suggestive Above 10 GeV, up-going /i's are used to determine z^-fluxes. The variation of 
1 up-going //-fluxes with the arrival direction can be used to study the oscillation parameters, since the distance to the 
place where v's are produced is determined by the zenith ang le of the arrival direction |9| |l0| [pd|. 

Atmospheric z/-fluxes have been calculated by Volkova ]l2f| , Mitsui et al. ]l3|| , Butkevich et al. |14| a nd Lipari |T^ 1 
mainly for high energies (from around 1 GeV to above 100,000 GeV). Gaisser et al. Jig , Barr et al. [17|], Bugaev and 
Naumov |l8| , Lee and Koh jl9), and Honda et al. ||(| calculated precisely the atmospheric i/-flux for low energies 
(< 3 GeV) . A calculation of low energy atmospheric z/-flux using the /x-flux observed at high altitudes has also been 
made §|. 

In this paper, we use essentially the same models for particle interaction, atmospheric structure, and cosmic ray 
fluxes with Ref. |2(J. The calculation method, however, is different for low energy (30 MeV - 3 GeV) and high energy 
(1 - 3,000 GeV) atmospheric u's. We employ a full Monte Carlo method at low energies, but use a hybrid method at 
high energies. 

The difficulties in the calculation of atmospheric (/-fluxes differ between high and low energies. In case of low 
energy z/s, the primary fluxes of cosmic ray components are relatively well known. However, the low energy cosmic 
ray fluxes (< 30 GeV) are modulated by solar activity, and are affected by the geomagnetic field through the rigidity 
(= momentum/charge) cutoff. For high energy z/s (> 100 GeV), the > 1,000 GeV cosmic ray flux is relevant. At 
these energies, solar activity and rigidity cutoff do not affect cosmic rays, but details of the cosmic ray flux are not as 
well-measured. Details of the hadronic interactions of cosmic rays with air nuclei are also a source of the uncertainty 
in the calculated i/-fluxes. At low energies, the proton-nucleus interaction at < 30 GeV is important. There have been 
many accelerator experiments studying hadronic interactions in this energy region, however, not many are suitable for 
our purpose. In the high energy proton-proton interactions, it is normally assumed that the spectrum of secondary 
particles satisfies the Feynman scaling hypothesis, which is confirmed by collider experiments up to a lab. energy of 
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3,000 TeV. Although there is a weak breaking of the hypothesis in the central region, it has no significant effect on 
atmospheric influxes. 

We employ an one-dimensional approximation in which all the secondary particles and the z/s keep the direction 
of their parent cosmic rays, throughout the energy range of concern. For high energy i/'s, this is a good approxi- 
mation because of the nature of hadronic interactions. We also expect that the directional-average of ^-fluxes may 
be calculated with good accuracy even at low energies. However, when we need information about the variation 
with direction of the low energy atmospheric i^-flux, especially for near horizontal directions, a three-dimensional 
calculation is necessary. At low energies, secondary particles are produced with large scattering angles by hadronic 
interactions, and the curvature of low energy /z's due to the geomagnetic field becomes sizable. We note that a three- 
dimensional calculation of ^-fluxes with the Monte Carlo method requires an enormous computation time, since we 
need to calculate the z^-flux at every position on the Earth, and for all direct ions with good statistics. 

Section |l| is devoted to the problems with low-energy primary cosmic rays (II A), such as the effect of solar activity 
(1IB) and the rigidity cutoff due to the geomagnetic field ( pi C| ). Also in [II D| , the primary cosmic ray fluxes are 
compiled for each chemical composition in the energy region of 100 GeV - 100 TeV, for use in the calculation 
of atmospheric z/s. The processes which take place during the propagation of co smic rays in the atmosphere are 
explained in section [II. The hadronic interaction model we employ is explained in [II A. The decay of mesons, such 
as 7r's and -ftf's, which are created in cosmic ray interactions, is the main source of the atmospheric z/s. These decay 



processes are summarized in section IIIB with a discussion of muon p olariz ations. In section [V, we ex plain the 



calculation of atmospheric ^-fl uxes, and the results are summarized in IV B for 30 MeV - 3 GeV, and in IV C for 
1 GeV - 3,000 GeV. In section [VD, atmosphericu-fluxes are calculated by the same method as the ^-fluxes, and 
are compared with the observed data. In section M, the uncertainties in the calculation of atmospheric influxes are 
discussed. In section VI , the major results of this work are summarized. 
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FIG. 1. 



Observed fluxes of cosmic ray protons, helium nuclei, and CNOs from the compilation of Webber and Lezniak [22]. 
(Dots represent data from Refs [23] and [24], diamonds from Ref. [25], crosses from Ref. [26], minuses from Ref. [27], 
squares from Ref. [28] , upward triangles from Ref. [29] , and downward triangles from Ref. [30] . ) Solid lines are our 
parametrization for solar mid., dash lines for solar min., and dotted lines for solar max. 
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II. THE FLUX OF COSMIC RAYS 



A. The Flux of Cosmic Rays below 100 GeV 

Primary cosmic ray fluxes are relatively well known in the low energy region (< 100 GeV), by which the low 
energy atmospheric ^-fluxes (< 3 GeV) are mainly created. However, the fluxes are affected by solar activity and 
the geomagnetic field. The effect of solar activity is known as the solar-modulation of cosmic rays, and is commonly 
parametrized by the sun-spot number or the count rate of neutron monitors. The effect of the geomagnetic field is 
represented as the rigidity cutoff of cosmic rays. In the following, the treatment of these effects in this calculation is 
explained. 

Webber and Lezniak have compiled the energy spectrum of primary cosmic rays for hydrogen, helium, and CNO 
nuclei in the energy range 10 MeV ~ 1,000 GeV p2| for three levels of solar activity (Fig. [j]). A similar compilation 
has been made by Seo et al. |3l) for hydrogen and helium nuclei, which agrees well with that of Webber and Lezniak. 
Seo et al. estimated that uncertainties in the instrumental efficiency (~ 12%) and exposure factor (2 - 3%) result in 
the overall uncertainty of the primary cosmic ray fluxes being ~ 15% ifBlfl . 

From the compilation of Webber and Lezniak, the chemical composition of cosmic rays is H(proton) ~ 90.6 %, He~ 
9.0 %, and CNO nuclei - 0.4 % above - 100 MeV/nucleus, and H - 95.2 %, He- 4.5 %, and CNO nuclei - 0.3 % 
above ~ 2 GeV/nucleus. The portion of other components (Ne, S, Fe, • • •) is so small that they can be neglected in 
the calculation of low-energy atmospheric ^-fluxes. It is noted that atmospheric v's are created through the hadronic 
interactions of cosmic rays and air nuclei, and therefore are dependent on the number of nucleons rather than the 
number of nuclei. The contribution of a heavier cosmic ray nucleus to the atmospheric j/-flux is larger than that of 
a cosmic ray proton. Helium nuclei carry ~ 15 % of the total nucleons in the cosmic ray flux and the CNO group 
carries ~ 3.6 % above ~ 2 GeV/nucleon. These effects are amplified by the effect of the geomagnetic field through the 
rigidity cutoff, since the rigidity for those nuclei is two times larger than protons with the same momcntum/nucleon. 
The details are given below. 



B. Solar Modulation 



The flux of low energy cosmic rays is modulated by solar activity. The solar wind drives back the low energy cosmic 
rays which are entering into the solar sphere of influence, and the strength of the solar wind varies with solar activity. 
This effect is more evident in the lower energy cosmic rays: the flux difference at solar maximum and solar minimum 
is more than a factor of two for 1 GeV cosmic rays, and it decreases to ~ 10 % for 10 GeV cosmic rays (Fig. Q). 

The primary flux for various levels of solar activity is parametrized by Nagashima et al. as a function of the 
count rate of a neutron monitor N by 

f(E k )dE k = ^u- 2 - 585 M(p, N)dE k , (2.1) 

where i stands for the kind of nucleus (=H, He, CNO ..), p for rigidity in GV (=GeV/c/Z), E k for kinetic energy per 
nucleon in GeV, and u for total energy per nucleon in GeV. The absolute flux value for each component is determined 
by 7i , where 7 H = 10.85 x 10 3 m 2 sec- 1 sr- 1 GeV _1 , 7h = 5.165 x 10 2 nAec^sr^GeV" 1 , and 7cno = 3.3 x 10~ 2 
m 2 sec _1 sr _1 GeV _1 respectively. The function M(p, N) is the modulation function defined by 

and V is the count rate of the neutron monitor at Mt. Washington with N max — 2465 count/hour. We take N — 2445 
for solar min., V = 2300 for solar mid., and N = 2115 for solar max. The results of the parametrization are shown 
in Fig. [I] by the solid, dashed and dotted lines, which agree well with data except for low energy helium nuclei 
(< 10 GeV/nucleon). However, this produces only a very small effect on the calculation of atmospheric z/'s due to the 
proton dominance of the cosmic ray flux and the relatively small contribution of cosmic rays of this energy. 



C. Rigidity Cutoff of the Geomagnetic field 

The geomagnetic field determines the minimum energy with which a cosmic ray can arrive at the Earth. This 
effect is caused by the magnetic shield effect for low energy cosmic rays. For the cosmic ray nucleus, the minimum 
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energy of cosmic rays arriving at the Earth is determined by the minimum rigidity (rigidity cutoff) rather than the 
minimum momentum. We note that the rigidity cutoff is a function of the entering position on the Earth and arrival 
direction (zenith angle, 8, and azimuth angle, <f). Since the mass/charge ratio of helium and CNO nuclei is twice that 
of protons, those nuclei carry lower energy nucleons into the atmosphere than protons. 

The actual geomagnetic field is represented by a multipole expansion of the spherical harmonic function as 
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with the potential function 
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Here R is the radius of the Earth and P™(x) is the associated Legendre function. The expansion coefficient (g™, h™) 
is compiled and reported by IAGA Division I working Group 1 p3| . As the geomagnetic field varies slowly with time, 
the coefficient is reported yearly with the time differential values. 

The value of the rigidity cutoff for the actual geomagnetic field can be obtained from a computer simulation of 
cosmic ray trajectories. In this simulation, an antiproton, which has the same mass as a proton but the opposite 
charge, is used as the test particle. We note that the change e <-> — e is equivalent to the change of t «-> — t in the 
equation of motion of a charged particle in a magnetic field; 

| = evx B, (2.5) 

where p is the momentum, v is the velocity, and B is the magnetic field. To determine the rigidity cutoff at different 
positions and for different directions, we launch antiprotons from the earth, varying the position and direction. When 
a test particle with a given momentum reaches a distance of 10 times the Earth's radius, where the strength of the 
geomagnetic field decreases to the same level as the interstellar magnetic field (~ 3 x 10~ 8 T), it is assumed that the 
test particle has escaped from the geomagnetic field. Assuming the momentum distribution of cosmic rays is isotropic 
in angular space, some cosmic rays which have the same rigidity can arrive at the Earth following the same trajectory 
but in the opposite direction. The rigidity cutoff is calculated as the minimum momentum with which the test particle 
escapes from the geomagnetic field. We note that for protons the rigidity and the momentum is the same quantity. 
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FIG. 2. 

The contour map of cutoff-rigidity for the v arrival directions at Kamioka. Azimuth angles of 0°, 90°, 180°, and 
270° show directions of south, east, north, and west respectively. 



In an one-dimensional approximation, we need the rigidity for the arrival direction of z/s at ^-detector sites. We 
found that the magnetic field calculated up to 5th order of the expansion gives almost the same result for rigidity cutoff 
as calculations with higher order expansions. The rigidity cutoff at Kamioka is shown as a contour map in Fig. 0. 
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We note that the dipole approximation is good near the magnetic equator. However, the multipole effect becomes 
important near the magnetic pole. For the calculation of atmospheric (/-fluxes, we need the multipole expression for 
the geomagnetic field even when the detectors are not located near the pole, since v's created near the pole arrive at 
the detector through the Earth. 

D. Cosmic Ray Fluxes above 100 GeV 

Cosmic rays with energy greater than 100 GeV, which are responsible for > 10 GeV atmospheric ^-fluxes, are not 
affected by solar activity or by geomagnetic effects. However, ther are few measurements of the cosmic ray chemical 
composition at these energies, especially above 1 TeV. Here, we compile the available data of cosmic ray fluxes for H, 
He, CNO, Ne-S, and Fe-group nuclei up to 100 TeV/nucleon. Above 100 TeV/nucleon, the cosmic ray spectrum is 
measured by the air shower technique and almost no direct measurements of cosmic ray particles are available. 

In Fig. ||, observed cosmic ray fluxes from Refs (2f| [Q |3(| |37j Q are summarized for H and He. We fitted the 
observed flux for > 100 GeV with a single power function: 

Flux(E) = A ■ {E/100 GeV) 7 , (2.6) 

and show the result in the same figure. We note that the data for He of Ryan et al. j25| are more than two times 
smaller than those of other groups and their error bars are larger than others. Therefore, we have not used their data 
in this analysis. The observed cosmic ray flux for CNO, Ne-S, and Fe-group nuclei from Refs [^6| |35| [B8| |3iJ pKj 
flOf [[ll| are shown in Fig. |], also with our fitted spectra. The parameters A, 7 are summarized in Tabic |for H, He, 
CNO, Ne-S, and Fe-group nuclei. 
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FIG. 3. 

The observed cosmic ray fluxes for H and He nuclei. The lines show our fits, as explained in the text. Diamonds 
are from Ref. [25], squares from Refs [34] and [35], upward triangles from Ref. [36], downward triangles for He from 
Ref. [37] , and crosses for H from Ref. [38] . The dashed line show the JACEE spectrum for H above 40 TeV [35] . 

For the calculation of the atmospheric ^-fluxes of energy region 1 - 3,000 GeV, we employed the superposition 
approximation. This approximation treats a nucleus as a bundle of independent nucleons, and considers the event 
caused by the nucleus as the sum of independent events caused by these nucleons. Therefore, we need the flux of each 
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nucleon rather than that of each chemical composition. In Fig. ||, we depicted the nucleon flux (proton and neutron) 
calculated from the single power function fit for H, He, CNO, Ne-S, and Fe-group cosmic rays. For < 100 GeV, 
we grouped the experimental points in Figs [j], [jj, and || in several energy bins and calculated the average flux over 
different groups. Thus calculated proton flux below 100 GeV agrees well with the parametrization of equation |2.l| 
with N = 2300 (solar mid.), while the neutron flux is larger. This difference has already appeared in Fig. |l|, but 
produces very small effects on the atmospheric z/-flux as explained before. It is noted that the proton cosmic rays 
constitute ~ 80 % of all nucleons at 100 GeV/nucleon, and helium ~ 15 %. This percentage decreases with energy, 
however, ~ 80 % of nucleons are carried by proton and helium cosmic rays even at 100 TeV/nucleon (Fig. ^|). The 
heavier nuclei are still a minor component of cosmic rays at this energy. 



FIG. 4. 



Observed cosmic ray fluxes for CNO, Ne-S, and Fe nuclei. The lines show the fitted result as is explained in the text. 
Pluses are from Ref. [26], minuses from Ref. [30], squares from Ref. [35], crosses from Ref. [38], dots from Ref. [39], 
upward triangles from Ref. [40], downward triangles from Ref. [41], and diamonds from Ref. [42]. 

Above 100 TeV, almost no observations of the cosmic ray chemical composition are available. However, the all- 
particle flux is well measured by the air-shower technique (e.g. see Ref. J43|). Converting the energy-per-nucleon 
spectra back to energy-per-particle, the H, He, CNO, Nc-S, and Fe-group nuclei fluxes were summed and compared 
with the observed all-particle spectrum |||] Q Q jSJ Q in Fig. 0. It is seen that the extrapolation of the calculated 
all-particle flux agrees well with the experimental data even at energies above 100 TeV. We used the nucleon flux 
calculated here up to 1,000 TeV (upper extrapolation for protons in Fig ||). The JACEE group suggested that the 
spectrum of cosmic protons becomes steeper (7 = 3.22) above 40 TeV [|4j (See Fig. ||). The effect of this steepening 
has been studied for {/-fluxes at 1,000 GeV, but is very small. Errors in the cosmic ray nucleon spectra obtained here 
are estimated to be < 10 % at around 100 GeV and increase to ~ 20 % at 100 TeV. 
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FIG. 5. 

The cosmic ray nucleon flux summed over all chemical components, and extrapolated (dashed line) up to 1,000 TeV. 
The lower line for protons above 40 TeV assumes the steepening suggested by the JACEE group. 
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Ratios of nucleons carried by H, He, CNO, Ne-S, and Fe-group nuclei. 
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FIG. 7. 



The all-particle flux summed over fitted spectra for H, He, CNO, Nc-S, and Fe-group nuclei compared with the 
observed data. Large squares from Ref. [35], downward-triangles from Ref. [36], small squares are from Ref. [43], 
minuses from Ref. [44], and upward-triangles from Ref. [45]. 



III. THE INTERACTION OF COSMIC RAYS IN AIR 



A. Hadronic Interaction 



As cosmic rays propagate in the atmosphere, they create 7r's and K's in interactions with air nuclei. These mesons 
create atmospheric i/'s when they decay as follows: 
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A cr + A a 



7T~ — > 



e+ + ^ e + P., 



A*" 



v.. 



(3.1) 



e + v e + v u 



The interactions of cosmic ray protons and nuclei with air nuclei are hadronic in character. We employed the 
NUCRIN Eg] M onte Carlo code for hadronic interactions for lab. energies < 5 GeV, and LUND code - FRITIOF 
version 1.6~|1] and JETSET version 6.3 El - for 5 GeV < E lah < 500 GeV. Above 500 GeV, the original code 
developed by Kasahara et al. (COSMOS) Q was used. The K/tt ratio is taken 7 % at 10 GeV, 11 % at 100 GeV, 
and 14 % at 1,000 GeV in lab. energy. We compared the output of those codes with available experimental data. 
Although there have been many accelerator experiments of proton-nucleus and Helium-nucleus interactions for < 
30 GeV, not many are applicable to our present purpose. For higher energies, the data of collider experiments are 
available for p — p and p — p collisions. We estimate the error of the atmospheric v calculation resulting from these 
Monte Carlo codes is around 10 %. 

Results from the LUND code were compared with the ir production spectrum in a cone of < 7.28° in p—Be collisions 
(24 GeV/c) [)5l| (Fig. 0). NUCRIN results were compared with the 7r production spectrum at 2.5° in p— C collisions 
(5 GeV/c) |52|| (Fig. The agreement of the LUND code and experimental data is quite good (~ 10 %) except for 
the very low momentum region (< 5 GeV). Since the energy spectrum of cosmic rays is steep, the spectrum of 7r's 
production by nucleons in the lower momentum region is relatively unimportant. The agreement of the NUCRIN code 
and experimental data seems not as good as that of the LUND code. The authors of the NUCRIN code estimated 
the disagreement of their output and experimental result as 10 - 20 % (47[| . In low energy cosmic ray interactions, the 
detailed structure of the 7r production spectrum may be smeared due to the flattening of the cosmic ray spectrum at 
low energies. 
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FIG. 8. 



Comparison of tt -production spectrum in p — Be collisions at p = 24 GeV/c between experiment (Ref. [50]) and 
the LUND Monte Carlo code. The direction of the 7r's are limited to < 7.28° from the direction of projectile protons. 
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FIG. 9. 

Comparison of Tr^-production spectrum in p— C collisions at p = 5.1GeV to the direction of 2.5° between experiment 
(Ref. [51])and the NUCRIN Monte Carlo code. P* denotes the momentum in the total center of mass system. 

For high energies, it is difficult to get experimental results of nucleus-nucleus or proton-nucleus interactions. How- 
ever, there are many results available from p-p(p) collider experiments. In the Fig. [lO], we present the experimentally 
determined pseudo-rapidity distribution for s/s = 53, 200, 546, and 900 GeV, and calculated results from the COS- 
MOS and LUND codes for the same energies and y/s = 30.6 GeV, corresponding to lab. energy of 500 GeV. Above 
V^s = 53 GeV, the results of the COSMOS code agrees with the experimental results within < 5 %. Also the agreement 
of the COSMOS and LUND codes is good at *Js = 30.6 GeV. However, the pseudo-rapidity distributions calculated 
by the LUND code is lower than the experiment results and those of the COSMOS code near rj = 0. Accordingly 
the multiplicity (the number of particles created by the interaction) in the LUND code is smaller than the experi- 
mental value above this energy. Therefore, above 500 GeV in lab. energy we used the COSMOS code for hadronic 
interactions. 




-4.0 -2.0 0. 2.0 4.0 

Pseudorapidity (T| ) 

FIG. 10. 



Pseudo-rapidity distributions calculated by the COSMOS code (solid line). Upper triangles, downward triangles, 
and squares represent the data from UA-5, SppS (Ref. [52]) for ^/s = 200, 546, and 900 GeV respectively. Crosses 
are data from UA-5, ISR (Ref. [53]) for y/s = 53 GeV. The dotted lines show the calculated results by LUND code 
for y/i = 30.6 and 53 GeV. 

In the COSMOS code, nucleus - nucleus interactions are treated as follows; first, the projectile (cosmic ray) nucleus 
fragments into smaller mass number nuclei and nucleons with given probabilities. Second, each fragmented nucleon 
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interacts with the target independently, and the sum of created particles in each interaction is considered as the 
product of the nucleus nucleus interaction. We take the average number of nucleons which interact with the target 
(air) nucleus as; 



< N >= A -^ (3- 2 ) 



where, A is the mass number of projectile nucleus, a™ el (<7™ ei ) is the inelastic cross-section of the proton (nucleus) 
with the air nucleus. We used the full treatment of COSMOS for the calculation of ^-fluxes in the 30 MeV - 3 GeV 
range, and the s uperposition model in the 1 - 3,000 GeV range. The validity of the superposition model is discussed 



later (Sec. |V C[) 



A qualitative discussion of the hadronic interactions is useful in order to understand the and v e excess over their 
antiparticles. In hadronic interactions of low energy where one tt production is dominant, the projectile charge is most 
often carried by the leading particle. Thus, the most probable interaction for cosmic ray protons of this energy is; 

P + A air -> n + 7T+ + X (3.3) 

where p stands for a proton and n for a neutron. In the higher energy region where protons cause the multiple 
production, there also is a similar effect that the 7r + 's produced in this interaction statistically have a larger energy 
than 7T~'s (see Fig. ^). Integrating with the steep spectrum of cosmic rays, there is an excess of 7r + 's over 7r~'s by ~ 
20 % for 7r-decay. Consequently, we expect ~ 20% excess of v e over D e . We note the neutron component of cosmic 
rays, which is carried by helium and other nuclei, creates a 7r _ 's excess over tt +, s from isospin-symmetry. However, 
this gives a small effect on the 7r-spectrum and the atmospheric neutrino flux due to the proton dominance in the 
cosmic ray flux. 



B. Decay of Mesons 

Neutrinos are mainly produced in the following decays of /i's, 7r's, and ICs [pq 
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(27 %) 


{K^ v ) 




-> 7r ± e =F J/ e (z? e ) 


(38.6 %) 





(3.4) 



The decay of charged 7r's and subsequent /i decay (tt — fi decay): 

i (3.5) 

/i ± —> e ± + v e (y) e + ^(fju), 

is dominant among these processes. Charmed particles, such as D and D also create ^'s, however, the contribution 
of charmed particle to atmospheric z/s becomes sizable only for E v > 100 TeV, which is far beyond the energy region 
of concern here. 

When a 7r ± decays at rest, the energy carried by v^ip^) is (m\ — m 2 ^)/2m 7T ~ 30 MeV, and the /x ± carries the 
rest of the energy. If we ignore the spin of the /i's, each decay particle, e^, u e {y e ) 1 and v^(i>^) carries 1/3 of fi^'s 
energy (~ 37 MeV) on average in the three body decay. When 7r's decay in flight, the tt energy is approximately 
divided into 1/4 to each decay product in the tt — [i decay on average. Thus, we expect the flux ratio to be roughly 
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(v e + v e )/{v^ + v^) = 1/2 and v^jv^ = 1 irrespective of the 7r spectrum. When the energy of the /x's becomes high 
(> 5 GeV), however, /i's tend not to decay in the air but to reach the Earth before decaying. In this case, /i's lose 
their energy in the Earth, and decay after they are almost stopped, or are captured by nuclei in the Earth. This 
effect reduces the ratio [y e + u e )/(v^ + 9^) above this energy. Also the ratio v^jv^ decreases with energy in the same 
energy region, since cosmic ray protons create more 7r + 's than 7r~'s, and there is a corresponding excess of /i + 's over 
/x~'s. We note that fj, polarization and /i energy loss in air are important for the precise calculation of atmospheric 
u's. Each has > 5 % effects on the u's energy created by /i's. The energy loss of /i's in the air is taken into account 
by the Monte Carlo method. The treatment of /i polarization is explained in the following. 

In the decay of charged it's, the resulting /i ± is fully polarized against (toward) the direction of fi motion in the 
charged 7r or K rest frame. In the subsequent /i decay, v e (v e ) is emitted to the forward direction of the /i's motion 
from the conservation of helicity. Thus the u e (u e ) resulting from /i decay has a larger energy than v^iyS) in the 7r rest 
frame. Since 7r's decay in flight, the polarization is not full. In general, the direction distribution of v e is proportional 
to (l + (cos8) where C is the polarization (£ — 1 is full polarization) and 8 is the angle between directions of ^ e -motion 
and /i-spin in the CM frame of the /i's. Practically, it is the polarization in the observer's frame that is important. 
Using the spin direction 3-vector £, the polarization along the /i momentum direction is 



? p_ _ EE* - 7ff mg 
\p\ P-P* 



c=c-^= ~ r (3.6) 



in the observer's frame [p6| . Here E is the /i energy, p the /i momentum , and is the Lorentz factor of the tt in 
observer's frame. E* and p* denote the /i energy and the momentum in the tt rest frame respectively. 

The above discussion can also be applied to /i's created in the K ± — > /i ± t' A1 (z7 M ) decay. For /i's resulting from the 
Kz\ v decay, the discussion in Ref. Q is applied. The small angle scattering of /i's in the atmosphere reduces the \x 
polarization. This depolarization effect is also evaluated by Hayakawa [[56| as of the order of 21MeV/wp, where v and 
p are velocity and momentum of /i's respectively. Therefore, the depolarization effect may be negligible for /i's which 
produce u's of energies > 100 MeV. If we ignore the effect of /i polarization, the calculated energy of u e decreases by 
~ 5 % on average, therefore, the flux of z/ e is estimated to be smaller by ~ 10 % at 500 MeV and ~ 15% at 3 GeV. 
This is an important effect for low energy i/s, which are observed in underground detectors |59f|. 



IV. THE FLUX OF ATMOSPHERIC NEUTRINOS 

A. ^—fluxes and Atmospheric Density Structure 

We note first that the density structure of the atmosphere is important, because it is the reason for a large zenith 
angle dependence of influxes. We take the US standard (60| for the density structure of the atmosphere. The chemical 
composition of atmosphere is approximated by nitrogen 78 %, oxygen 21 %, and argon 1 % in our calculation. For 7r's 
and ICs propagating in the atmosphere, the decay and interaction with air nuclei are competitive processes. When 
the relation 

El , , 

CT 4.1 

mcr an 

is satisfied, both processes work at nearly the same rate. Here, r is the life time of the meson, E is the energy, a the 
interaction cross-section of the meson and air nuclei, and n the number density of air nuclei. This condition (4.1) is 
rewritten for the energy of mesons with which decay and interaction take place almost equally as 

2 f 12 (GeV, for 7r±) ] 
E ~ — = I 22 (GeV, for K9) } x P[sea level] , (4.2) 



CTon 



90 (GeV, for K d 



where P[ sea i eV ei] — 1.225kgm -3 . Most mesons which have smaller energies than that given by this equation decay. 
Slant entering cosmic rays interact with atmospheric nuclei at a higher altitude than vertically entering cosmic rays 
on average. Therefore, the decay probability is larger for the mesons created by slant entering cosmic rays than for 
those of vertical cosmic rays. We expect a larger i/-flux from near horizontal directions than from the vertical. 

The first interaction of vertical cosmic rays takes place at an altitude of 15 ~ 20 km, where the density of air is ~ 10 
times less than that at the sea level, since the interaction mean-free-path (MFP) for cosmic ray protons is ~ 100 g/cm 2 
in column density. Mesons with energy < 12 GeV, which create z/s of < 3 GeV, decay before interacting with air 
nuclei. However, for mesons with energy > 100 GeV, the probability of interaction becomes sizable. Therefore, the 
variation of atmospheric i/-fluxes with the zenith angle increases with the ^-energy. 
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B. Low Energy v Flux (30 MeV - 3 GeV) 



At low energies, although the zenith angle dependence of ^-fluxes caused by the atmospheric structure is not large, 
a significant directional variation is caused by the rigidity cutoff. In the one-dimensional approximation we adopted, 
we expect larger ^-fluxes from the low rigidity cutoff direction and a smaller ^-fluxes from the high rigidity cutoff 
direction. There should be a large directional dependent variation of ^-fluxes for Kamioka (Fig. ||) , especially a large 
deficit from east-horizontal direction. In the actual case, however, it may be difficult to observe those variations. 
The direction of a v is different from that of the parent cosmic ray, because the mesons are produced with slightly 
different directions from that of incident particle, and again the i/'s are created with slightly different directions from 
those of the mesons in their decay. We expect this effect to be rather small, but there is another smearing effect of 
direction in the ^-detector. When a low energy v (< 3 GeV) creates a lepton by a quasi-elastic process, the lepton 
has a typical angle of 50 — 60° from the v direction Jjl). Therefore, it may be difficult to observe the deficit of ^-flux 
from east-horizontal direction in Kamioka. Thus the direction dependence of atmospheric z/-flux is small for lower 
energy ;/s, especially when they are observed in the detector. We present here the atmospheric i^-flux averaging over 
all directions. 

For low energy ^-fluxes, we employed a full Monte Carlo method. The calculation itself is r ather straightforward; 
first, the nucleus and the primary energy of the cosmic ray are sampled with equation (2.1). Second, the arrival 
direction is sampled uniformly. When the rigidity is smaller than the rigidity cutoff, the cosmic ray is discarded. 
When the rigidity is larger than the rigidity cutoff, the cosmic ray is put to the propagation code of cosmic rays in 
the atmosphere, which controls the interaction of cosmic rays, the decay of secondary particles, and the energy losses 
in the atmosphere. The COSMOS code controls all these processes. 

The results are summarized in Table |l| for Kamioka and in Table || for the 1MB site, both for solar mid. Also 
in Fig. [[ll i/ft + and v e + v e fluxes are depicted for Kamioka at solar mid., and compared with the other results. 
Since ^-fluxes for solar mid. are not available for other authors, we averaged solar max. and solar min. values. Flux 
differences between solar max. and solar min. are ~ 8 % at 100 MeV and ~ 3 % at 1 GeV for Kamioka, and ~ 
12 % at 100 MeV and ~ 4 % at 1 GeV for 1MB due to its lower cutoff-rigidity. We note that these calculated results 
are smoothly connected to the atmospheric neutrino fluxes calculated by the hybrid method for high energies (1 - 
3,000 GeV) at around 3 GeV, where the effect of the rigidity cutoff is small. 
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The atmospheric z/-fluxes multiplied by E^, for the Kamioka site at solar mid. (solid lines). BGS are from Ref. [16], 
BN from Ref. [17], and LK from Ref. [18]. The dotted line is the result from the calculation for high energy without 
the rigidity cutoff, and averaged over all directions. For details, see the next section (IV C). 

In Fig. |l2|, we show the flux ratio by ^-species along with those of other authors. We note that although the 
calculation method and some of physical assumptions are different among these authors, the ratio (y e + D e )/(v^ + v^) 
is very similar each other. The relatively large difference in D e /v e among them may reflect the difference of calculation 
scheme and/or the physical assumptions. 
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FIG. 12. 



The flux ratio of ^-species calculated for Kamioka. BGS are from Ref. [16], BN from Ref. [17], and LK from 
Ref. [18]. as before. 



C. High Energy v Flux (1 - 3,000 GeV) 

For the calculation of atmospheric z/s at high energies (1 - 3,000 GeV), we employed the superposition model 
and a hybrid method. We note that for helium nuclei, the complete treatment of COSMOS and the superposition 
model give essentially the same result. Let us measure the atmospheric depth by the column density from the top of 



the atmosphere. With equation (3.2), the average number of nucleons which interact with an air nucleus in a depth 
interval [x, x + dx] is given by 

< iV > e x a — = A- -e x a =A-— --e x a — (4.3) 
Xa <ja Xa X p Xa 

per nucleus. The nucleus MFP in the column density of air is denoted by Xa, and the proton MFP denoted by X p . 
We note that a oc 1/A. Since the fragmentation of helium nuclei into deuterons occur with only a small probability, 
we can safely assume that helium nuclei always fragment into four nucleons {p,p, n, n}. The non-interacting nucleons 
at the first collision also interact with air nuclei in succeeding processes. The average number of nucleons which did 
not interact in the helium collision at the depth x\ and that do interact in [x, x + dx] is given by 

f x i i * dx 

(A- < N >)[ —e x a ■ — e a p dxi] dx = A ■ [e x ? ~ e x a] — (4.4) 
Jo Xa X p X p 



for a nucleus. Adding (4.3) and (4.4), we get the average number of interacting nucleons in a depth interval [x, x + dx] 



as 

- -2- dx . . - * - » . dx . - dx ... 
A ■ e x a — + A ■ [e x p - e x a ] — = Ae x p —, (4.5) 

A p Ap Ap 

which is exactly the same with the case of A-independent nucleons. We note that up to E cr = 100 TeV, around 80 % 
of all nucleons are carried by hydrogen and helium nuclei. 



13 



For heavier cosmic rays than helium, the probability of fragmentation into smaller nuclei (e.g., helium nucleus = 
a particle, etc.) is large at collisions. Therefor e, th e average number of nucleons which interact with the target (air) 



nucleus is smaller than that given by equation (|3.2|). Thus, the interaction height (depth) distribution of each nucleon 
is different from proton cosmic rays in general. This distribution is represented by the effective MFP of nucleon for 
heavier nuclei, and is calculated by the COSMOS code as ~ 100 g/cm 2 for CNO and ~ 94 g/cm 2 for Fe nuclei. The 
MFP is 87.4 g/cm 2 for protons, 36.6 g/cm 2 for CNO's, and 15.6 g/cm 2 for Fe nuclei at 1 TeV/nucleon. Thus the 
effective MFP for nucleons carried by those nuclei is more similar to the proton MFP than to those nuclei. These 
facts reasonably justify the superposition model, even for nuclei heavier than heliums. 

In our hybrid method, we calculate the z^-yield function for protons, T^,{E, E v , 0), which denotes the number of v 
in the energy region from E v to E v + dE u , created by a proton with energy of E incident from zenith angle 9 with 
the Monte Carlo method. We executed the Monte Carlo simulation for proton cosmic rays with energies from 1 GeV 
to 1,000 TeV in Alog(£ p ) = 0.2 steps, and for zenith angle bins of cosfl = 0.0 - 0.1, 0.1 - 0.2, .... 0.9 - 1.0. For each 
zenith angle bin, 300,000 protons were injected to the Monte Carlo code at 1 GeV, 100,000 protons at 10 GeV, 
1,000 protons at 100 TeV, and 300 protons at 1,000 TeV. 

The proton z/-yield function for v^'s and zVs calculated by the Monte Carlo method is fitted by a function: 

with parameters {A,B,C}. This fit agrees with the result of the Monte Carlo method very well in the region 



3 ?S E p /E u < 10 3 , as is seen in Fig. O. The solid lines show our fit with equation 4.6. As the cosmic ray spectrum 
is approximately proportional to -E^ in a wide energy range, the quantity shown in Fig. [l^ roughly stands for the 
relative contribution to a fixed energy z/s from the cosmic ray in a logarithmic energy bin. The largest contribution 
to fixed energy atmospheric i/'s comes from the cosmic ray with E cr ~ lOE^ both for i^'s and u e 's. We note that the 
contribution from E v jE v < 3 and E p /E v > 10 3 is very small(< 5 %). The yield functions of u e 's and P e 's decreases 
more rapidly than u^s and C^'s, and they change their shape for E v > 100 GeV. 
For the neutron incident, we assumed 

r£(E,E v ,6)=ri»(E,E v ,0) and 

r£(E,E v ,6)=Ti3(E,E v ,8) . (4.7) 



This assumption is justified for z/s produced through 7r's, but not for those produced through A"'s. The K~ / K + ratio 
is rather an universal quantity for p-A and A-A interactins (~ 0.8) at high energies S|. However, we note that the 
portion of proton nucleon is still ~ 80 % in the cosmic ray at 100 TeV (Fig. [j]). Therefore, the assumption leads to 
a maximum of 10 % errors in the K~ / K + -r&tio. We note that this assumption affects only on the vjv ratio at high 
energies. We expect almost no effect for vjv ratio at low energies, because of the proton dominance in the low energy 
cosmic ray. Also the v + z/-flux is not affected in any energy region. 

The atmospheric z^-fluxes were calculated by integrating those v— yield functions with the nucleon fluxes shown in 
Fig | 

pEmax 

F„{E U ) = / [r£{E,E v ,ff) x F P (E) +rg {E,E V , 6) x F n (E)]dE (4.8) 

J Emin 

and 

pEmax 

F^{E„)= [r£(E,E u ,9)xF p (E) + rfi(E,E u ,6)xF n (E)]dE , (4.9) 

J Emin 

where F v (Fv) is the atmospheric i/(i>)-flux, F p and F n are the proton and neutron-fluxes respectively, and v stands 
for or v e . We took E min equal to the z/-energy and E max — 1, 000 TeV. 



14 



V H -Vertical 



v e -Vertical 




E P /E V 



FIG. 13. 



E P /E v 



^-yield function for protons with fixed E v multiplied by E 17 for (left) and v e (right) for near vertical directions 
(cos# = 0.9 — 1.0). Monte Carlo results are shown by squares, and fitted results by lines. 



The results of the integrations (4 
calculated fluxes below 



and (4.9) are summarized in Tables 
10 GeV smoothly connect to the low energy calcul. 
shows t hat th e systematic difference between the calculations for high energy and low energy atmospheric ^-fluxes 
(section IV B) is very small. With the limited number of protons used in the calculation of v-y\e\& function at high 



fv|, [V]jv| and |vn] down to 1 GeV. The 
ation at around 3 GeV (Fig. |ll|). This 



energies, the errors involved in {^4, B, C} become large due to large fluctuations. Therefore, the error on the flux 
value increases at high energies. This error is estimated to be < 5 % for and below 100 GeV, and increases 
to 10 % at around 3, 000 GeV for near horizontal directions. The error for v e and v e fluxes is larger than and 
fluxes; it may exceed 10 % at 1,000 GeV for near horizontal directions and at 100 GeV for near vertical directions. 
For v^s and C^'s, however, we extrapolate the flux value up to 3.16 x 10 4 GeV. We note that large errors at these 
energies do not affect much the expected flux of up-going /i's, because the contribution of v's with E v > 1,000 GeV 
is estimated less than 15 %. The value in the parentheses is estimated to have a larger error than 10 %. 

The ^-fluxes below 3 GeV could differ substantially from the true flux, since solar modulation and geomagnetic 
effects have not been included. The fluxes for < 3 GeV in these tables should be regarded as those for solar mid. 
and low cutoff rigidity (< 3 GV). They also are depicted for the near vertical and near horizontal directions in FigJIJ 
with the results of other authors. The v^iv^ fluxes are almost proportional to E~ z from 1 GeV to 1 TeV. On the 
other hand, the v e {ve) fluxes decrease, being proportional to -fiT 3 ' 5 or steeper above 10 GeV. Our ^-flux for near 
horizontal directions agrees with that of Volkova |0| and Lipari |Te| above 10 GeV to within ~ 5 %. For near vertical 
directions, our calculation is 10 - 15 % larger than that of Volkova and Mitsui et al. |l3), and by 5 - 10% than that 
of Lipari. For i/ e -fluxes, our calculation is larger than others by 10-20% above 10 GeV for both directions. 

If we take the lower line for the extrapolation of the proton flux in Fig. || for > 40 TeV, the and z? M -fluxes 
decrease ~ 10 % at 3,000 GeV, but the y de crease onl y a few % at 1,000 GeV. Even if we reduce the upper end of 
the integrations (E max ) to 100 TeV in (4.8) and (4.S), i.e. if there was a sharp cutoff of the cosmic ray spectrum 
at 1,000 TeV, the and P M fluxes are reduced by only ~10 % at 1,000 GeV. Therefore, even if the cosmic ray 
composition above 100 TeV is quite different from our assumption, our calculation does not give a different result for 
the atmospheric v— fluxes below 1,000 GeV by more than 10 %. However, for the accurate calculation of atmospheric 
i/s above 1,000 GeV, it is necessary to determine the flux and composition of primary cosmic rays above 1,000 TeV 
per nucleon accurately. 
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To study the portions of v's resulting from if's and from 7r— u-decay, we also calculated the ^-fluxes produced only 
with 7T — /z-decay and showed the ratio to the total flux in Fig. |l5|. The contribution of charmed mesons to atmospheric 
z/s is very small in this energy range. The /z-flux which resulted from 7r-decay was also calculated and the ratio to 
the total flux is also shown in the same figure. It is seen that the i/'s created by tt — fj, decay are the minor component 
above ~ 30 GeV for near vertical directions and ~ 500 GeV for near horizontal directions both for v^'s and v e 's. It 
is also seen that although the main source of //-fluxes is 7r-decay up to 10,000 GeV, the main source of atmospheric 
i/'s is -FT's above ~ 30 GeV for near vertical directions and above ~ 500 GeV for near horizontal directions. 
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v n + and v e + v e fluxes for near horizontal (cos# = — 0.1) and near vertical (cos# = 0.9 — 1.0) directions. 
Volkova is from Ref. [12], Mitsui et al from Ref. [13], Butkevich et al from Ref. [14], and Lipari from Ref. [15]. We 
used the fluxes of cos 9 = 0.05 for near horizontal and cos 9 = 1 for near vertical directions of Volkova and Lipari, 
those of 9 zen ith = 87° for near horizontal and 9 zen i t h — for near vertical directions of Mitsui et al, and those of 
cos 9 = for near horizontal and cos 9 = 1 for near vertical directions of Butkevich ct al. 



The ratios v^jv^ and v e jv e are shown in in Figs ^ and ^ for near vertical and near horizontal directions with 
results of other authors. It is difficult to determine the ratio accurately with the Monte Carlo method, especially near 
t he h igh energy end. We do not consider the variation of ~ 0.05 for each ratio to be meaningful. Also the assumption 
(4.7) is the source of errors in the v^jv^ and v e jv e ratios at high energies, since the main source of atmospheric 
v's is the -ftT-decay. However, the i^/i^-ratio is clearly smaller than other calculations |l3] fl5|| , especially for near 
vertical directions even in the energy region where the the assumption (4.7) is good (< 30 GeV). For the v e /v e ~ 
ratio, the agreement with other calculations is better than the i/^/P^-ratio, but also shows significant differences at 
E v < 30 GeV for near vertical directions. Since we have sufficient statistics at E v < 30 GeV both for and v e in 
our Monte Carlo method, we consider that those differences result from differences in the calculation schemes, such 
as the interaction m odel a nd/or the atmospheric density structure. Our ratios calculated here agree well with those 
calculated in Section IV B at ~ 3 GeV. 
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The ratio of + u^, v e + v e , and ^+ + fluxes produced by 7r — /i-decay to the total fluxes. Solid lines show 
that for + t^, dashed lines for v e + D e , and dotted lines for ^+ + [i~ . Near horizontal denotes the average over 
cos 9 = — 0.1 and near vertical denotes that over cos6> = 0.9 — 1. 
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This work 



Butkevich et al 
Lipari 



10 u 10 10 10' 

Neutrino Energy (GeV) 



z/^/z^-ratio. Notations are the same as Fig. [14]. 
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FIG. 17. 



^ e /z/ e -ratio. Notations are the same with Fig. [14]. 

The ratio {y^ + v^)l{y e + D e ) calculated here is also compared with other authors in Fig. [T|. Our result is smaller 
than others for > 10 GeV both for near vertical and near horizontal directions, reflecting the larger v e -&uxes than 
other calculations in this energy region. This may be due to differences in the calculation scheme. However, all results 
show good agreement in the < 10 GeV region, except for that of Mitsui et al. for near horizontal directions. Their 
result is larger than others by ~ 50 % in the < 10 GeV for near horizontal directions, and by 10 — 15 % even for near 
vertical directions in the same energy region. We note that Mitsui et al. did not take into account the effect of muon 
polarization. This explains the difference for near vertical directions, but not for near horizontal directions. There 
seems to be other differences in their calculation scheme. 




FIG. 18. 



(lip + v^) J '(y e + j/ e )-ratio. Notations are the same with Fig. [14]. 



D. Flux of Atmospheric Muons 



There have been many observations of /i-flux for near vertical and near horizontal directions at sea level. Since the 
/^-flux is the complementary part of the ^-flux in ir — ji decay, the calculation of ^-fluxes is often examined by the 
comparison of the calculated /i-flux with the observed flux. We calculated the atmospheric /^-flux using the hybrid 



method; we first calculate the /i-yield function, then integrated it with the nucleoli flux. In Fig. 19, we present our 
calculated /i + + fj," fluxes for near vertical directions (cos 9 — 0.9 — 1) and the observed near vertical flux, and in 
Fig. [2^ for near horizontal directions. We note that although the main source is different mesons for atmospheric /i's 



18 



and i/s at high energies, the nucleon energy dependence of the /i-yield function is very similar to that of the z^-fiux 
above 100 GeV for all zenith angles. 

In the same figures, we show the calculated results of Butkevich et al. [[l4j and Lipari We used the result of 
cos 9 = 1 for the near vertical direction and of cos 9 = for near horizontal directions for these authors. Since the 
DEIS group measured the ^-flux at many zenith angles, we summed their data in the zenith angle bins of 78° — 84° 
and 84° - 90° corresponding to cos 9 = 0.1 - 0.2 and cos 9 = - 0.1. The MUTRON group measured the muon flux 
in the direction of 86° — 90° and the flux averaged zenith angle is 89° . 

It can be seen that the agreement of our calculation and the observations is good, although there is some variation 
in observed fluxes among the different groups. Our calculation agrees with DEIS data within ~ 15 % at all energies, 
and < 5 % for 100 - 1,000 GeV. The agreement is especially good for the near horizontal directions. Taking into 
account the fact that average zenith angle of the MUTRON data is 89°, the agreement of MUTRON data and our 
calculation is also very good. However, since the main parent mesons are different for z/s and /z's, these agreements 
shown above do not fully justify our calculation of ^-fluxes. 
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FIG. 19. 



Calculated /i + + /i~~flux for cos 9 = 0.9 — 1 and observed fluxes. Dots are from Ref. [62], squares from Ref. [63], 
crosses from Ref. [64], minuses from Ref. [65], upward triangles from Ref. [66], downward triangles from Ref. [67], 
vertical diamonds from Ref. [68] , and horizontal diamonds from Ref. [69] . Also shown are the calculated results from 
Ref. [14] (dashed thin line) and Ref. [15] (dotted thin line) both for cos9 = 1. 
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Calculated /i + + /j -flux for cos6> = — 0.1 (solid line) and 0.1 — 0.2 (dotted line), and the observed flux by the 
DEIS (Ref. [69]) and MUTRON (Ref. [70]) groups. Also shown are calculated results from Ref. [14] (dashed thin 
line) and Ref. [15] (dotted thin line) both for cos 9 = 0, for comparison with the MUTRON data. We note that the 
DEIS observation is added for the zenith angle bins of 78° — 84° and 84° — 90° corresponding to cos# = 0.1 — 0.2 and 
cos 6 = 0-0.1. 



The charge ratio /i + /fi~ was also calculated and shown in Fig. ^l] with the result of other authors |1J] As 
we do not consider that the variation of < 0.05 to be meaningful in our calculation, a constant value of ~ 1.25 can 
explain our results both for near vertical and near horizontal directions. Othe r ca lculations show an increase with 
energy, especially for the near vertical directions. We note that the assumption (4.7) is almost valid for /i + //z~-ratio, 
since the main source is still the 7r — /i-decay, and the experimental results do not show such an increase but are 
consistent with the constant value ~ 1.25 for both directions. 
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FIG. 21. 

Comparison of calculated /i + //Lt~-ratios with observations. For the near vertical directions, dots are from Ref. [62], 
crosses from Ref. [69] , squares from Ref. [72] , minuses from Ref. [73] , upward triangles from Ref. [74] , and diamonds 
from Ref. [75]. For the near horizontal directions, triangles are Ref. [70] and squares from Ref. [71]. Also for near 
horizontal direction, calculated /i + //i~-ratios are averaged for cos# = — 0.2 for the results of Lipari and this work. 



At low energies, the production altitude of v's is 80 - 300 g/cm 2 for vertical directions, which corresponds to 9 - 
18 km (Fig. ^2|) altitude. Since /i's travel ~ 6 km before decay at 1 GeV on average, fi's observed at sea level are not 
directly related to v's at low energies (< 1 GeV). For the examination of the calculated atmospheric ^-flux at low 
energies, the observation of /^-fluxes at the production height is necessary. 
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Production height distribution of u's with energy > 400 MeV by vertical cosmic rays. 



Recently, the /i~-fiux has been measured with good accuracy by the MASS (Matter-Antimatter Superconducting 
Spectrometer) experiment |7(J at high altitudes. In Fig. we compared the observed //"-fluxes and our calculation 
at the same altitudes. This //-flux was calculated by the full Monte Carlo method, the same as the low energy v-ftux. 
The agreement of experiment and calculation is very good except for very high altitudes (< 37 g/cm 2 ) and the low 
momentum region at low altitudes (615 g/cm 2 ). We note that since the low energy atmospheric i/'s are created at 
the altitude of 80 - 300 g/cm 2 , the contribution of atmospheric z/s created at very high altitude (< 37 g/cm 2 ) or low 
altitude (615 g/cm 2 ) to the total z^-flux is relatively small. (See Fig. p2|.) 
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FIG. 23. 



A comparison of \x -fluxes observed by the MASS (Matter-Antimatter Superconducting Spectrometer) experi- 
ment [76] at several altitudes and this work. 



V. THE SYSTEMATIC ERROR AND OTHER UNCERTAINTIES 



The systematic error in the atmospheric ^-fluxes comes mainly from the uncertainty of the cosmic ray primary 
flux. Even at low energies, where the primary cosmic ray flux is rather well studied, it is difficult to determine the 
absolute value due to the uncertainties in the instrumental efficiency (~ 12%) and exposure factor (2 - 3%) 
These uncertainties in the primary cosmic ray flux increase with energy. In our compilation, the error in the fit is 
~ 10 % for the nucleon flux at 100 GeV and ~ 20 % at 100 TeV. Assuming ~ 10% uncertainty below 100 GeV, the 
systematic error in the atmospheric ^-fluxes is estimated to be ~ 10% at < 3 GeV, increasing to ~ 20 % at 100 GeV, 
and remain almost constant up to 1,000 GeV. 

We note that the uncertainty of the primary cosmic ray flux increases more rapidly than the fitting error above 
10 TeV/nucleon. The JACEE group suggests a steepening of the cosmic ray proton spectrum above 40 TeV pi| . 
Using this steep proton spectrum, the atmospheric i/-ffux decreases 3-4 % at 1,000 GeV, ~ 10 % at 3,000 GeV, and 
> 20 % above 10,000 GeV. Below 300 GeV, the difference is negligibly small. 

The uncertainty in the nucleon spectrum above 100 TeV, where almost no direct observations are available, is 
crucial for the calculation of atmospheric j/'s above 1,000 GeV, The air shower technique, which is commonly used 
to study cosmic rays above 100 TeV, can not determine the chemical composition with the accuracy we need. The 
nucleon spectrum could be very different from that assumed here, depending on the major chemical component of 
cosmic rays above 100 TeV. 

The interaction model is another source of systematic errors. In our comparison, the agreement of the LUND model 
and the COSMOS code with the experimental data is < 10 %. The agreement of the NUCRIN code and experimental 
data is not as good as the LUND code. The authors of the NUCRIN code claim that the agreement is within 10 - 
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20 % [J47|]. However, the hadronic interaction below 5 GeV is related to very low energy atmospheric ^-fluxes; when 
we switch off the hadronic interaction for < 5 GeV, they decrease 2 - 5 % at 1 GeV, 15 - 25 % at 300 MeV, and 45 - 
55 % at 100 MeV dependeng on the rigidity cutoff. Therefore, we may conclude that the error of the NUCRIN code 
does not affect the atmospheric i/-flux by more than 10 % in our calculation for > 300 MeV. The systematic error 
caused by the hadronic interaction model is estimated to be ~ 10 % above 300 MeV. 

The calculation scheme - the one-dimensional approximation in all energy regions and the superposition model 
for nucleus-nucleus interactions at high energies - is also the possible source of some systematic error. Since one- 
dimensional approximation is justified at high energies, it is expected to be accurate above 3 GeV. With the one- 
dimensional approximation, however, the calculation of rigidity cutoff is very simplified, and this may result in a 
systematic error in the absolute value of the atmospheric i^-fluxes of 10 - 20 % at 100 MeV and 5 % at 1 GeV. Other 
effects caused by the one-dimensional approximation and errors due to the superposition model at high energies are 
considered to be small compared to other errors. 

The calculation method also results in an error in the atmospheric j/-flux. With the full Monte Carlo calculation, 
it is estimated to be < 2 — 3 % for 30 MeV - 3 GeV due to the statistics. With the hybrid method, the statistics and 
the fitting error are sources of the uncertainty. Both the statistics and fitting error are combined and estimated to be 

< 5 % up to 100 - 300 GeV for and D^, and up to 30 - 100 GeV for v e and v e , depending on the zenith angle. 
Combining all the systematic and non-systematic errors, the total error is estimated as 20 % at 100 MeV, 15 % 

from 1 GeV to 100 GeV, and 20 - 25 % at the highest energy in our calculation. However, the errors of the species 
ratio is smaller than the absolute value, since the j'-species ratio is not affected much by the uncertainty of primary 
fluxes and the calculation scheme. It is estimated to be < 10 % below 100 GeV for v/v and < 5 % below 30 GeV for 
(Up + v ll )/(v e + D e ). These errors also increase to 10 - 15 % at the highest energies in our calculation. 

Although the main parent meson is different for high energy atmospheric /i's and i/'s, one may consider that the 
comparison of calculated and observed atmospheric //-flux reduces the systematic error due to the primary cosmic 
ray flux. The agreement of our calculation, DEIS data, and MUTRON data for near horizontal directions (Fig. |20| ) 
suggests that the systematic error of the atmospheric ^-flux caused by the uncertainty of the cosmic ray flux may be 

< 10 % at E v — 100 - 1,000 GeV in our calculation. However, we note that there are similar problems for observation 
of /i-fluxes with that of primary cosmic rays: determining the instrumental efficiency and the exposure factor. The 
MASS group claims that their total uncertainty in efficiency is around 20 %. This systematic error for ground based 
experiments could be smaller, but it is seen in Fig. |l9| that //-fluxes observed for near vertical directions by different 
groups differ by more than 20 %. Also taking the possible systematic error in the iif/7r-ratio into account, we conclude 
that the systematic error caused by the uncertainty in the primary flux is ~ 20 % at E u = 1, 000 GeV. 



VI. SUMMARY 



In this paper, we presented the calculation of atmospheric ^-fluxes as follows; first, we summarized the physics 
related to the primary cosmic ray flux. The low energy cosmic ray flux was parametrized following the work of 
Nagashima et al. |52| based on the compilation of the cosmic ray spectrum by Webber and Lezniak |22| for solar max., 
mid., and min. In order to calculate the rigidity cutoff due to the geomagnetic field, we simulated the trajectories of 
cosmic rays. Also the primary cosmic ray flux in the 100 GeV - 100 TeV range was compiled for H, He, CNO, Ne-S, 
Fe nuclei groups calculating the nucleon flux for protons and neutrons. We used the NUCRIN interaction model ^] 
for < 5 GeV, LUND model (JETSET ver. 6.3) @ || for 5 - 500 GeV, and COSMOS for > 500 GeV for 
the hadronic interaction interactions of cosmic rays. The atmospheric ^-flux was calculated with a full Monte Carlo 
method for 30 MeV - 3 GeV, and with a hybrid method for 1 GeV - 1,000 GeV. 

One of the most important results for low-energy ^-fluxes is that the ratio (v e + v e )/(v^ + v^) is almost the same 
(~ 0.5) as other calculations for 100 MeV - 3 GeV jlTj |l8| |ll| pi]], whereas underground detectors found a significant 
difference in the ratio (e-like event)/ (//-like event) from the value expected from the calculated atmospheric ^-flux. 
We note that the quantity {y e + v e )l{v^ + v^) remains relatively unaffected by variations in the interaction model 
and primary cosmic ray spectrum including the chemical composition. The difference between the observed and the 
expected value of the ratio (e-like event)/(/x-like event) might be explained by other physics, such as ^-oscillations. 

If this difference is to be explained by ^-oscillations with Am 2 ~ 10~ 2 eV 2 , the up-going //-flux, which is induced by 
the high energy ^-flux, will show a different zenith angle dependence from the expected one. We calculated the zenith 
angle dependence of the atmospheric ^-flux in detail: for each zenith angle bin of cos 8 = — 0.1,0.1 — 0.2, • • • 0.9 — 1.0, 
for 1 - 1,000 GeV. This atmospheric ^-flux could be used to calculate the expectation flux of up-going /t's. We note 
that the absolute value and the ratio are connected to the lower values smoothly at ~ 3 GeV. Also they are compared 
with the calculation of other authors (l2) |]| @ |]| . 
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Atmospheric /z-fluxes were also calculated at sea level and at high altitudes. They were compared with the experi- 
mental data and the agreements are found to be satisfactory. The agreement at high altitude is especially important 
for the calculation of the atmospheric i/-flux at low energies. Although the parent mesons are different for atmo- 
spheric j/s and yti's, our calculation, DEIS data, and MUTRON data agree very well each other for the /z-flux of near 
horizontal directions at high energies (> 100 GeV). We conclude that we have used a reasonable primary cosmic ray 
spectrum, chemical composition, and interaction model. 

We stress again that the main source of the systematic error in the atmospheric f-flux is the uncertainty of the 
primary cosmic ray flux. Especially for the calculation of the atmospheric i/-flux above 1,000 GeV, the lack of 
knowledge of the cosmic ray flux above 100 TeV/nucleon is crucial. 
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TABLE I. Fitted parameters for each chemical component of cosmic rays. 
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TABLE III. The low energy u-Rux for 1MB (solar mid., m" 2 sec" 1 sr" 1 GeV" 1 ) 
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TABLE IV. i^— flux xE'l (m 2 sec 1 sr 1 GeV 2 ) calculated with the hybrid method. The value above 1 x 10 4 is the smooth 



extension with a power fit. 
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344.6 


314.5 


293.7 


276.6 


261.9 


249.3 


238.7 


251.2 


525.2 


444.3 


384.9 


341.6 


310.8 


289.0 


271.2 


256.2 


243.5 


233.0 


316.2 


521.3 


440.6 


380.9 


337.1 


305.5 


282.9 


264.5 


249.0 


236.2 


225.8 


398.1 


516.4 


435.1 


374.9 


330.6 


298.4 


275.2 


256.3 


240.7 


227.8 


217.4 


501.2 


510.3 


427.2 


366.1 


321.4 


289.0 


265.7 


246.8 


231.3 


218.6 


208.3 


631.0 


503.0 


417.3 


355.0 


309.9 


277.6 


254.6 


236.1 


220.9 


208.5 


198.5 


794.3 


494.6 


405.4 


341.7 


296.4 


264.5 


242.3 


224.5 


209.8 


197.8 


188.2 


1000. 


485.3 


391.7 


326.5 


281.1 


250.0 


228.9 


212.0 


198.1 


186.7 


177.4 


1259. 


474.4 


375.7 


308.8 


263.5 


233.4 


213.7 


198.1 


185.1 


174.3 


165.3 


1585. 


461.6 


356.9 


288.2 


243.2 


214.3 


196.5 


182.2 


170.2 


160.0 


151.5 


1995. 


447.1 


335.8 


265.4 


220.9 


193.7 


177.8 


165.0 


154.1 


144.6 


136.5 


2512. 


431.1 


313.0 


241.2 


197.6 


172.2 


158.3 


147.1 


137.3 


128.6 


120.9 


3162. 


406.0 


296.5 


228.2 


185.2 


158.4 


142.0 


128.8 


(117.5) 


(107.8) 


(99.44) 


10 4 


(301.2) 


(188.1) 


(127.1) 


(93.0) 


(73.6) 


( 62.5) 


(53.6) 


(46.1) 


(39.6) 


(34.1) 


3.16 • 10 4 


(201.2) 


(99.5) 


(55.45) 


(34.8) 


(24.7) 


(19.4) 


(15.4) 


(12.1) 


(9.5) 


(7.4) 
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TABLE V. ^-flux xE'l (m 2 soc x sr 1 GeV 2 ) calculated with the hybrid method. The value above 1 x 10 4 is the smooth 



extension with a power fit. 



_E„(GeV)\cos6> 


.0 - .1 


.1 - .2 


.2 - .3 


.3 - .4 


.4 - .5 


.5 - .6 


.6 - .7 


.7 - .8 


.8 - .9 


.9 - 1. 


1.000 


217.1 


210.4 


203.3 


196.0 


188.3 


180.6 


173.5 


167.0 


161.1 


155.7 


1.259 


246.6 


237.4 


228.0 


218.4 


208.7 


199.1 


190.4 


182.6 


175.7 


169.6 


1.585 


274.4 


262.2 


249.9 


237.7 


225.6 


213.8 


203.4 


194.5 


186.9 


180.5 


1.995 


299.6 


283.5 


267.9 


252.8 


238.2 


224.3 


212.4 


202.6 


194.6 


188.2 


2.512 


323.6 


302.6 


283.2 


265.0 


248.2 


232.7 


219.7 


209.2 


200.9 


194.6 


3.162 


348.0 


321.3 


297.5 


276.5 


257.7 


241.1 


227.5 


216.4 


207.7 


201.1 


3.981 


371.7 


338.6 


310.5 


286.6 


266.3 


249.1 


235.0 


223.6 


214.5 


207.5 


5.012 


393.4 


353.9 


321.5 


295.0 


273.4 


255.9 


241.6 


229.9 


220.5 


213.2 


6.310 


412.6 


366.6 


330.2 


301.4 


278.9 


261.5 


247.1 


235.3 


225.6 


218.0 


7.943 


429.2 


377.0 


336.7 


305.9 


282.7 


265.4 


251.1 


239.2 


229.4 


221.4 


10.00 


443.5 


385.1 


341.2 


308.4 


284.5 


267.3 


253.2 


241.2 


231.2 


223.0 


12.59 


455.6 


391.0 


343.5 


308.9 


284.3 


267.3 


253.3 


241.3 


231.2 


222.7 


15.85 


465.5 


394.8 


343.7 


307.4 


282.3 


265.6 


251.7 


239.7 


229.4 


220.5 


19.95 


473.2 


396.2 


341.8 


303.9 


278.4 


262.1 


248.5 


236.5 


225.9 


216.6 


25.12 


478.4 


395.3 


337.8 


298.5 


272.8 


257.0 


243.6 


231.6 


220.8 


211.1 


31.62 


481.3 


393.0 


332.8 


292.5 


266.7 


251.3 


238.2 


226.3 


215.3 


205.2 


39.81 


482.1 


389.9 


328.0 


286.9 


261.2 


246.1 


233.3 


221.4 


210.3 


200.0 


50.12 


480.9 


386.3 


323.3 


281.9 


256.2 


241.5 


228.9 


217.1 


206.0 


195.6 


63.10 


477.6 


382.0 


318.7 


277.3 


251.8 


237.3 


224.9 


213.2 


202.2 


191.8 


79.43 


472.2 


377.1 


314.2 


273.2 


247.9 


233.5 


221.3 


209.8 


198.9 


188.7 


100.0 


464.8 


371.6 


309.8 


269.4 


244.4 


230.2 


218.1 


206.8 


196.2 


186.2 


125.9 


456.0 


365.2 


304.8 


265.2 


240.6 


226.5 


214.5 


203.4 


193.0 


183.4 


158.5 


446.2 


357.5 


298.4 


259.6 


235.4 


221.5 


209.7 


198.8 


188.7 


179.3 


199.5 


435.4 


348.6 


290.8 


252.8 


229.1 


215.4 


203.8 


193.1 


183.2 


174.0 


251.2 


423.8 


338.6 


282.0 


244.8 


221.6 


208.2 


196.9 


186.4 


176.6 


167.6 


316.2 


411.5 


327.7 


272.1 


235.7 


213.2 


200.1 


189.1 


178.8 


169.2 


160.3 


398.1 


398.7 


316.4 


262.0 


226.4 


204.3 


191.5 


180.6 


170.4 


160.9 


152.0 


501.2 


385.6 


305.5 


252.4 


217.5 


195.7 


182.7 


171.6 


161.3 


151.6 


142.6 


631.0 


372.3 


294.9 


243.2 


209.1 


187.2 


173.8 


162.3 


151.6 


141.7 


132.5 


794.3 


358.9 


284.6 


234.6 


201.0 


179.0 


164.9 


152.8 


141.6 


131.3 


121.9 


1000. 


345.5 


274.7 


226.4 


193.2 


170.9 


156.0 


143.1 


131.3 


120.7 


110.9 


1259. 


331.6 


264.6 


217.9 


185.2 


162.4 


146.5 


132.7 


120.4 


109.3 


(99.4) 


1585. 


317.1 


253.8 


208.7 


176.3 


153.0 


136.0 


121.4 


108.5 


(97.2) 


(87.1) 


1995. 


302.1 


242.4 


198.8 


166.6 


142.8 


124.8 


109.5 


(96.2) 


(84.7) 


(74.7) 


2512. 


286.7 


230.6 


188.2 


156.3 


132.0 


113.2 


(97.4) 


(83.9) 


(72.4) 


(62.6) 


3162. 


278.89 


210.88 


165.26 


134.27 


113.09 


(98.3) 


(85.8) 


(74.9) 


(65.5) 


(57.2) 


10 4 


(206.9) 


(133.8) 


(92.1) 


(67.4) 


(52.6) 


(43.2) 


(35.7) 


(29.4) 


(24.1) 


(19.6) 


3.16 • 10 4 


(138.2) 


(70.8) 


(40.2) 


(25.3) 


(17.6) 


(13.4) 


(10.2) 


(7.7) 


(5.8) 


(4.3) 
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TABLE VI. ^ e -flux xE't (nr 2 scc _1 sr _1 GoV 2 ) calculated with the hybrid method. 



£„(GeV)\cos0 


.0- .1 


.1 - .2 


.2 - .3 


.3 - .4 


.4- .5 


.5 - .6 


.6 - .7 


.7- .8 


.8 - .9 


.9 - 1. 


1.000 


131.2 


124.3 


116.7 


109.0 


101.3 


94.0 


87.2 


81.0 


75.6 


70.9 


1.259 


150.9 


140.8 


130.5 


120.4 


110.7 


101.7 


93.6 


86.4 


80.2 


74.9 


1.585 


166.3 


154.8 


142.2 


129.5 


117.2 


105.9 


95.9 


87.4 


80.4 


75.0 


1.995 


179.2 


165.4 


150.1 


134.7 


120.1 


106.9 


95.5 


86.0 


78.6 


73.1 


2.512 


192.5 


172.0 


152.6 


134.9 


119.1 


105.5 


94.0 


84.6 


77.1 


71.3 


3.162 


205.4 


176.7 


152.7 


132.7 


116.2 


102.5 


91.2 


81.9 


74.4 


68.4 


3.981 


216.8 


181.7 


153.6 


131.1 


113.0 


98.4 


86.6 


77.1 


69.4 


63.2 


5.012 


226.2 


186.9 


155.3 


129.9 


109.7 


93.5 


80.7 


70.5 


62.5 


56.3 


6.310 


233.6 


190.4 


155.5 


127.5 


105.3 


87.8 


74.2 


63.7 


55.6 


49.7 


7.943 


239.0 


190.2 


151.9 


122.1 


99.1 


81.4 


67.9 


57.6 


50.0 


44.3 


10.00 


242.3 


186.2 


144.9 


114.2 


91.5 


74.5 


61.8 


52.3 


45.3 


40.1 


12.59 


243.5 


178.8 


134.8 


104.3 


82.8 


67.3 


55.9 


47.6 


41.4 


36.8 


15.85 


242.5 


168.3 


122.5 


93.0 


73.4 


59.9 


50.4 


43.4 


38.2 


34.1 


19.95 


239.1 


156.4 


109.8 


82.0 


64.6 


53.1 


45.3 


39.7 


35.4 


31.8 


25.12 


233.3 


144.7 


98.4 


72.4 


56.9 


47.2 


40.8 


36.2 


32.7 


29.5 


31.62 


225.2 


133.1 


88.0 


63.9 


50.2 


42.0 


36.7 


33.0 


30.1 


27.3 


39.81 


215.2 


121.9 


78.6 


56.5 


44.4 


37.4 


33.1 


30.1 


27.6 


25.1 


50.12 


203.6 


111.2 


70.2 


50.1 


39.5 


33.5 


29.9 


27.5 


25.3 


22.9 


63.10 


190.7 


100.9 


62.7 


44.5 


35.1 


30.0 


27.1 


25.0 


23.1 


20.7 


79.43 


176.8 


91.3 


56.0 


39.6 


31.4 


27.0 


24.5 


22.8 


21.0 


18.7 


100.0 


162.3 


82.2 


50.0 


35.3 


28.1 


24.3 


22.2 


20.7 


19.0 


16.7 


125.9 


147.4 


73.6 


44.6 


31.5 


25.2 


22.0 


20.2 


18.8 


17.2 


14.9 


158.5 


132.6 


65.7 


39.7 


28.2 


22.7 


19.9 


18.3 


17.1 


15.5 


13.1 


199.5 


118.1 


58.4 


35.4 


25.3 


20.5 


18.1 


16.7 


15.5 


13.9 


(11.5) 


251.2 


104.1 


51.7 


31.6 


22.8 


18.6 


16.5 


15.2 


14.0 


12.4 


(10.0) 


ft 


on q 
yu.y 


40.0 


08 1 
Zo. 1 


on k 
zu.o 


1 ft Q 


1 C 1 

10.1 


lo.y 


10 7 
1Z. / 


l~\ i n\ 
(ll.UJ 


{°-<) 


398.1 


78.5 


39.9 


25.0 


18.5 


15.4 


13.8 


12.7 


(11.5) 


(9.8) 


(7.5) 


501.2 


67.1 


34.9 


22.3 


16.7 


14.1 


(12.7) 


(11.6) 


(10.4) 


(8.7) 


(6.4) 


631.0 


56.8 


30.3 


19.8 


15.1 


12.9 


(11.6) 


(10.6) 


(9.4) 


(7.6) 


(5.4) 


794.3 


47.6 


26.2 


17.6 


13.7 


(11.8) 


(10.7) 


(9.8) 


(8.5) 


(6.7) 


(4.6) 


1000. 


39.5 


22.6 


15.7 


12.5 


(10.9) 


(9.9) 


(9.0) 


(7.7) 


(5.9) 


(3.8) 
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TABLE VII. P e -flux xEl (nr 2 scc _1 sr _1 GeV 2 ) calculated with the hybrid method. 



£„(GeV)\cos0 


.0- .1 


.1 - .2 


.2 - .3 


.3 - .4 


.4 - .5 


.5 - .6 


.6 - .7 


.7- .8 


.8 - .9 


.9 - 1. 


1.000 


109.2 


103.1 


96.5 


89.8 


83.2 


77.0 


71.1 


65.9 


61.3 


57.3 


1.259 


126.1 


117.3 


108.3 


99.6 


91.4 


83.7 


76.8 


70.7 


65.4 


61.0 


1.585 


139.2 


129.2 


118.4 


107.5 


97.1 


87.5 


79.1 


71.9 


66.1 


61.5 


1.995 


149.9 


138.1 


125.1 


112.0 


99.7 


88.7 


79.1 


71.2 


65.0 


60.5 


2.512 


160.7 


143.4 


127.1 


112.2 


99.1 


87.8 


78.3 


70.5 


64.3 


59.6 


3.162 


170.7 


146.7 


126.8 


110.3 


96.7 


85.4 


76.2 


68.7 


62.6 


57.8 


3.981 


179.2 


150.2 


127.2 


108.8 


94.1 


82.2 


72.7 


65.0 


58.8 


53.9 


5.012 


186.0 


154.0 


128.3 


107.8 


91.3 


78.2 


67.9 


59.7 


53.3 


48.3 


6.310 


191.2 


156.4 


128.2 


105.7 


87.7 


73.6 


62.5 


54.0 


47.5 


42.7 


7.943 


194.8 


155.8 


125.1 


101.2 


82.6 


68.3 


57.3 


48.9 


42.7 


38.2 


10.00 


196.8 


152.3 


119.2 


94.6 


76.2 


62.5 


52.1 


44.4 


38.6 


34.4 


12.59 


197.1 


146.0 


111.0 


86.5 


69.0 


56.4 


47.2 


40.3 


35.2 


31.3 


15.85 


195.8 


137.3 


100.8 


77.1 


61.3 


50.3 


42.4 


36.6 


32.2 


28.8 


19.95 


192.6 


127.6 


90.4 


68.1 


53.9 


44.5 


38.0 


33.3 


29.6 


26.5 


25.12 


187.6 


118.0 


81.1 


60.1 


47.5 


39.5 


34.1 


30.2 


27.0 


24.2 


31.62 


180.9 


108.6 


72.6 


53.2 


42.0 


35.1 


30.6 


27.3 


24.6 


21.9 


39.81 


172.6 


99.5 


65.0 


47.1 


37.1 


31.2 


27.4 


24.7 


22.3 


19.7 


50.12 


163.0 


90.6 


58.0 


41.7 


32.9 


27.8 


24.6 


22.2 


20.1 


17.7 


63.10 


152.3 


82.1 


51.6 


36.9 


29.1 


24.7 


22.0 


20.0 


18.1 


15.7 


79.43 


140.8 


74.0 


45.9 


32.6 


25.8 


22.1 


19.8 


18.0 


16.2 


13.9 


100.0 


128.8 


66.3 


40.7 


28.9 


22.9 


19.7 


17.7 


16.1 


14.4 


12.3 


125.9 


116.5 


59.0 


36.0 


25.5 


20.3 


17.6 


15.9 


14.5 


12.8 


(10.7) 


158.5 


104.4 


52.3 


31.8 


22.6 


18.1 


15.7 


14.2 


12.9 


11.4 


(9.4) 


199.5 


92.5 


46.1 


28.0 


20.0 


16.1 


14.0 


12.7 


11.5 


(10.1) 


(8.1) 


251.2 


81.1 


40.3 


24.6 


17.7 


14.3 


12.5 


11.4 


(10.3) 


(8.9) 


(7.0) 


Q1 (\ 9 


/ U.O 


oo. i 


zi.o 


10. D 


1Z. f 


119 
11. Z 


/in i)\ 
(.1U.2J 


(9.1) 


(7 8 s ! 


(b.U) 


398.1 


60.3 


30.4 


18.9 


13.8 


11.3 


(10.0) 


(9.1) 


(8.1) 


(6.8) 


(5.1) 


501.2 


51.2 


26.2 


16.5 


12.2 


(10.1) 


(9.0) 


(8.1) 


(7.2) 


(5.9) 


(4.3) 


631.0 


43.0 


22.4 


14.4 


(10.8) 


(9.0) 


(8.0) 


(7.3) 


(6.4) 


(5.1) 


(3-6) 


794.3 


35.8 


19.1 


12.5 


(9.5) 


(8.0) 


(7.2) 


(6.5) 


(5.6) 


(4.4) 


(3.0 


1000. 


29.4 


16.2 


(10.8) 


(8.4) 


(7.2) 


(6.4) 


(5.8) 


(4.9) 


(3.8) 


(2.5) 
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